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ABSTRACT 

Recent observations of thermally emitting isolated neutron stars revealed spectral fea- 
tures that could be interpreted as radiative transitions of He in a magnetized neutron 
star atmosphere. We present Hartree-Fock calculations of the polarization-dependent 
photoionization cross sections of the He atom in strong magnetic fields ranging from 
10^^ G to 10^^ G. Convenient fitting formulae for the cross sections are given as 
well as related oscillator strengths for various bound-bound transitions. The effects 
of finite nucleus mass on the radiative absorption cross sections are examined using 
perturbation theory. 
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1 INTRODUCTION 



An important advance in neutron star astrophysics in the 
last few years has been the detection and detailed studies 
I of surface emission from a large number of isolated neutron 
. stars (NSs), inclu ding radio pulsars , magnetars, and radio - 
■ quiet NSs (e.g., iKaspi et al.l l2006l : iHarding fc Lail 120061 ). 
, This was made possible by X-ray telescopes such as Chan- 
' dra and XMM-Newton. Such studies can potentially provide 
invaluable information on the physical properties and evolu- 
, tion of NSs (e.g., equation of state at super-nuclear densities, 
cooling history, surface magnetic field and composition). Of 
great interest are the radio-quiet, thermally emitting NSs 
fe.g.. lHaberil|2006[ ): they share the common property that 
their spectra appear to be entirely thermal, indicating that 
the emission arises directly from the NS surfaces, uncon- 
taminated by magnetospheric emission. The true nature of 
these sources, however, is unclear at present: they could be 
young cooling NSs, or NSs kept hot by accretion from the 
ISM, or magnetar descendants. While some of these NSs 
(e.g., RX J1856.5— 3754) have featureless X- ray spectrum 
rema rkably well described by blackbody fe.g.. iBurwitz et"ail 
I2OO3I ) or by emissi on from a cond ensed surface covered by 
a thin atmosphere ||Ho et al.ll2007l ). a single or multiple ab- 
sorption features at E ~ 0.2-1 keV have been dete cted from 
several sources (see Ivan Kerkwiik fc KaplanllioOTi ): e.g., IE 
1207.4-5209 (0.7 a nd 1.4 keV, possibly also 2.1, 2 8 keV ; 
ISanwal et all I2OO2I : iDe Luca et"aLl |2004|: iMori et al.ll2005l ). 
RX J1308.6-h2127 (0.2-0. 3 keV: iHaberl et al.l 12005). RX 
J1605.3-h3249 (0.45 keV : Ivan Kerkwiik et aiT 120041 ) . RX 
J0720.4-3125 (0.27 k eV: iHaberl et al.ll2006l ). and possibly 
RES 1774 (~ 0.7 keV: IZane et al.ll2005l ). The identifications 



of these features, however, remain uncertain, with sugges- 
tions ranging from proton cyclotron lines to atomic transi- 
tions of H, He. or m i d-Z atoms in a strong magnetic field (see 



Sanwal et all I2OO2I: IHo fc Lail [20oi : iPavlov fc Bezchastnovl 
20051 : iMori fc Hoi |2007| ). Clearly, understanding these ab- 



sorption lines is very important as it would lead to direct 
measurement of the NS surface magnetic fields and composi- 
tions, shedding light on the nature of these objects. Multiple 
lines also have the potential of constraining the mass-radius 
relation of NSs (through measurement of gravitational red- 
shift). 

Since the thermal radiation from a NS is mediated by 
its atmosphere (if T is sufficiently hig h so that the surface 
does not condense into a solid; s ee, e.g., Ivan Adelsberg et al.l 
I2OO5I : iMedin fc Lail |2006| . |2007| ). detailed modelling of ra- 
diative transfer in magnetized NS atmospheres is impor- 
tant. The atmosphere composition of the NS is unknown 
a priori. Because of the efficient gravitational separation 
of light and heavy elements, a pure H atmosphere is ex- 
pected even if a small amount of fallback or accretion oc- 
curs after NS formation. A pure He atmosphere results if 
H is completely burnt out, and a heavy-element (e.g., Fe) 
atmosphere may be possible if no fallback/accretion oc- 
curs. The atmosphere composition may also be afi'ected by 
slow) diffusive nuclear burnin g in the outer NS envelope 
Chang. Arras fc Bildsten|[2003 ). as well as by the bombard- 
ment on the surface by fast part icles from NS magneto- 
spheres (e.g., iBeloborodov fc Thompson 2007 ). Fully ion- 
ized atmosphere models in various magnetic field regimes 
have been extensively studied (e.g., IShibanov et al.l Il992l : 
IZane et allboOll : IHo fc Laill200lt). including the effect of vac- 
uum polarization (see lHo fc'Laill2003l : lLai fc Hdl20o3 . l2003l : 
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Ivan Adelsberg fc Laill2006l ). Because a strong magnetic field 
greatly increases the binding energies of a toms, mo lecules, 
and other bound species (for a review, see these 
bound states may have a ppreciable abundances in the NS 
atmos phere, as guessed bv lCohen. Lodenguai. fc Ruderman 



Thurner et al] Il993l: IJones et all Il999l : iMori fc Hailevl 



20021 : lAl-Huiai fc Schmelchej l2003bl l calculated oscillator 



197C ) and confirmed by calculations of 



19971 ') and IPotekhin. Chabrier fc Shibanov 



Lai fc Salpetei 
(|l999l ). Eaih 



considerations of partially i onized and strongly magn etized 
atmospheres fe.g.. iRaiagopal. Romani fc Milleiill997h relied 
on oversimplified treatments of atomic physics and plasma 
thermodynamics (ionization equilibrium, equation of state, 
and nonideal plasma efi'ects). Recently, a thermodynami- 
cally consistent equation of state and opacities for mag- 
netized [B = 10^^ — 10^^ G), partially ionized H p lasma 
have been obtained l|Potekhin fc Chabriedl2003l . |2004) . and 
the effect of bound atoms on the dielectric tensor of the 
plasma has also been studied ([Potckhin ct al. 200^. These 
improvements have been incorporate d into partially ion- 
ized, magnetic NS atmosph ere models (|Ho et al.ll2003l . [20071 : 
iPotekhin et"aLll2004 [2009 ) . Mid-Z element atm ospheres for 
10^^ - 10" G were recently studied by iMori fc Hoi 
|2003). 

In this paper we focus on He atoms and their radiative 
transitions in magnetic NS atmospheres. It is well known 
that for B ^ Z^Bq, where Z is the charge number 
of the nucleus and Bo = e^ml/Ti^c = 2.35 x lO'' G, 
the binding energy of an atom is significantly increased 
over its zero-field value. In this strong-field regime the 
electrons are confined to the ground Landau level, and 
one may apply the adiabatic approximation, in which 
electron motions along and across the field are assumed 
to be decoupled from each other (see Sect. l2.l| l. Using 
this approximation in combination with the Hartree-Fock 
method ("IDHF approximation"), several groups calcu- 
lated binding energies for the helium atom dProschel et al.l 



and molecules (INeuhauser, Langanke fc KooninI 19861: 


Neuhauser. Koonin fc Langanke 


1 19871: MiUer fc Neuhauser 


199l|: iLai, Salpeter fc Shapiro 


1992f). Mori fc Hailev 



l|2002l ) developed a "multiconfigurational perturbative 
hybrid Hartree-Fock" approach, which is a perturba- 
tive improvement of the IDHF method. Other methods 
of calculation include Th omas-Fermi-like models (e.g., 
[Abrahams fc Shapiro' 19911 ), t he d e nsity functional the - 

de^ 



ory (e.g., iRelovskv fc Rude: 
variational 



19961 : iMedin fc Lail |2006|) 



method s~ 7e.g- . iMiillerT Il984l: IVincke fc Bayel 



, . . - 

1 19891 : iJones et al.l 1 19991 : iTurbiner fc Guevaral |2006|) 
and 2D Hartree-Fock mesh calculations (|lvanovllT994 : 
llvanov fc Schmelcheij 12000*) which do not directly employ 
the adiabatic approximation. 

In strong magnetic fields, the finite nuclear mass and 
centre-of-mass mo tion affect the atomic structure in a non- 
trivial way (e.g., iLail I2OOII : see Sect.©. The stronger B 
is, the more important the effects of finite nuclear mass 
are. Apart from the H atom, these effects have been cal- 
culated only for the He atom which rests as a whole, but 
has a moving nucleus (Al-Huiai fc Schmelcher' '2003a"b'), 
and for the He"*" ion ([B ezchastnov. Pavlov fc Ventura 1998; 
IPavlov fc BezchastnovlFoOSi ). 

There were relatively few publications devoted to ra- 
diative transitions of non-hydro genic atoms in strong mag- 



netic fields. Several authors (|Miller fc Neuhauser Il99ll : 



strengths for bound-bound transitions; I Miller fc Neuhauser] 
(I99I) present ed also a few in t egrate d bound-free oscilla- 
tor strengths. iRaiagopal et al.l l| 19971 ) calculated opacities 
of strongly magnetized iron, using photoionization cross 
sections obtained by M. C. Miller (unpublished). To the 
best of our knowledge, there were no published calculations 
of polarization-dependent photoionization cross sections for 
the He atom in the strong-field regime, as well as the calcu- 
lations of the atomic motion effect on the photoabsorption 
coefficients for He in this regime. Moreover, the subtle effect 
of exchange interaction involving free electrons and the pos- 
sible role of two-electron transitions (see Sect. l3.2| l were not 
discussed before. 

In this paper we perform detailed calculations of radia- 
tive transitions of the He atom using the IDHF approx- 
imation. The total error introduced into our calculations 
by the use of these two approximations, the Hartree-Fock 
method and the adiabatic approximation, is of order 1% 
or less, as can be seen by the following considerations: The 
Hartree-Fock method is approximate because electron corre- 
lations are neglected. Due to their mutual repulsion, any pair 
of electrons tend to be more distant from each other than 
the Hartree-Fock wave function would indicate. In zero-field, 
this correlation effect is especially pronounced for the spin- 
singlet states of electrons for which the spatial wave function 
is symmetrical. In strong magnetic fields {B ^ Bq), the elec- 
tron spins (in the ground state) are all aligned antiparallel 
to the magnetic field, and the multielectron spatial wave 
function is antisymmetric with respect to the interchange of 
two electrons. Thus the error in the Hartree-Fock approach 
is expected to be less than the 1% ac curacy characteristic of 
zero-field Hartree-Fock c alculations (|Neuhauser et aj]|l987l: 
Schmelcher. Ivanov fc Be ckon 199^; for B = Q see Scrinzil 
1998^). The adiabatic approximation is also very accurate at 
B > Z^Bo- Indeed, a comparision of the ground-sta te en- 
ergy values calculated here to those of llvanovl l|l994 ) (who 
did not use the adiabatic approximation) shows an agree- 
ment to within 1% for B = 10^^ G and to within 0.1% for 
B = 10" G. 

The paper is organized as follows. Section [2] describes 
our calculations of the bound states and continuum states 
of the He atom, and section |3] contains relevant equations 
for radiative transitions. We present our numerical results 
and fitting formulae in section |3] and examine the effects of 
finite nucleus mass on the photoabsorption cross sections in 
section [S] 



2 BOUND STATES AND SINGLY-IONIZED 
STATES OF HELIUM ATOMS IN STRONG 
MAGNETIC FIELDS 

2.1 Bound states of the helium atom 

To define the notation, we briefly describe IDHF calcula- 
tions for He atoms in strong magnetic fields. Each electron 
in the atom is described by a one-electron wave function 
(orbital). If the magnetic field is sufficiently strong (e.g., 
B 2> 10^" G for He ground state), the motion of an electron 
perpendicular to the magnetic field lines is mainly governed 
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by the Lorentz force, which is, on the average, stronger than 
the Coulomb force. In this case, the adiabatic approximation 
can be employed - i.e., the wave function can be separated 
into a transverse (perpendicular to the external magnetic 
field) component and a longitudinal (along the magnetic 
field) component: 



(1) 



Here Wm is the ground-state Landau wave function (e.g.. 
iLandau fc Lifs hitz 1977) given by 



VF™(rx) = 



po\/27rm! \\/2po 



exp 



-P 



(2) 



where (p, kp) are the polar coordinates of rx , po = 
(Ticj eHY^^ is the magnetic length and jmv is the longitudi- 
nal wave function which can be calculated numerically. The 
quantum number m (^ for the considered ground Landau 
state) specifies the negative of the z-projection of the elec- 
tron orbital angular momentum. We restrict our considera- 
tion to electrons in the ground Landau level; for these elec- 
trons, m specifies also the (transverse) distance of the guid- 
ing centre of the electron from the ion, — (2m -f l)^''^po. 
The quantum number v specifies the number of nodes in 
the longitudinal wave function. The spins of the electrons 
are taken to be aligned anti-parallel with the magnetic field, 
and so do not enter into any of our equations. In addition, 
we assume that the ion is completely stationary (the 'infinite 
ion mass' approximation) . In general, the latter assumption 
is not necessary for t he applicability of the adiabatic ap- 
proximation (see, e.g.. |Potekhin|[l994l ). The accuracy of the 
infinite ion mass approximation will be discussed in Sect.jS] 
Note that we use non-relativistic quantum mechanics 
in our calculations, even when RoJse ^ m^c? ot B > Br 



Bo/a^ = 4.414 x lO" G. This is valid for two reasons: (i) 
The free-electron energy in relativistic theory is 



3Q 



E 



2 2 I 2 4 
C + m^c 



B 

bT, 



1/2 



(3) 



For electrons in the ground Landau level (wl =0), Eq. (O 
reduces to iJ ~ meC^ +p'i/{2me) for p^c <^ meC?\ the elec- 
tron remains non-relativistic in the z direction as long as 
the electron en ergy is much le ss than meC?; (ii) it is well 
known (e.g.. iSokolov fc Ternov|[l986i 'l that Eq. ^ describes 
the transverse motion of an electron with tll = at any field 
strength, and thus Eq. ((2)| is valid in the relativistic theory. 
Our calculations assume that the longitudinal motion of the 
electron is non-relativistic. This is valid for helium at all field 
strengths considered in this paper. Thus relativistic correc- 
tions to our calculated electron wave functions, binding ener- 
gies, and transition cross sectio ns are all small. Our ap prox- 
imation is justified in part by Chen fc GoldmanI l| 19921 ). who 
find that the relativistic corrections to the binding energy of 
the hydrogen atom are of order AS/S ~ 10"^ '^ - lO""'^ for 
the range of field strengths we are considering in this work 
{B = 10^2 - 10^* G). 

A bound state of the He atom, in which 
one electron occupies the (miJ^i) orbital, and the 
other occupies the (77121^2) orbital, is denoted by 
|mii^i, 77121^2) = I W^mi /mi 1^1, W^ma/ma 1^2) (clearly, 
I mi 1^1, 77121^2) ~ \m2V2 , rniUi)) . The two-electron wave 



function is 



[W/mi(rix)/mii,i(zi) 



XWm2 (T'2±)/m2 
-Wm2{r\l_)fm2V2{zi)Wmi{r21_)fmivi{z2)\ . (4) 

The one-electron wave functions are found using 
Hartree-Fock theory, by varying the total energy with re- 
spect to the wave functions. Th e total energy is given by 
(see, e.g.. iNeuhauser et al.lll987l ): 



E = Ek + EeZ + -Bdir + -Boxc , 

where 

-Ze^Yl j d2|/„.(2)|V™ 



Edir 



Ec 



E 



dzdz' |/mi.(s)|^ |/mv(2:')r 



4 E 



dzdz' /mV'(2)/m>'(z) 



and 

V;„(2) = y dr 



!IVm(rx)|' 



dr I dr 



dr I dr'i 



, \Wm{r 



\r' — r\ 

X 1^:., (rx)IVm(rx)W^:.(rl)W^w (r-1) 
Variation of Eq. ([Sjl with respect to fmu (z) yields 

Ze Vm[z) 



(5) 

(6) 
(7) 

(8) 

(9) 

(10) 
(11) 

(12) 



frnf {z) 



2me Az'^ 

+ e^ ^ J dz' \fm'v'{.z')fDram'{.Z - z) - £. 
ra ' v' 

= ^ y Az f*nu (z')/m',y {z')E,nm' {z " z')f^,^, {z) . 

in' v' 

(13) 

In these equations, asterisks denote complex conjugates, and 
fmi'(z) = Afmv/Az. The wave functions fmv{z) must satisfy 
appropriate boundary conditions, i.e., fmv ^ as a ^ ±00, 
and must have the required symmetry [fmv[z) = ±/m,^(— z)] 
and the required number of nodes {y). The equations are 
solved iteratively until self-consistency is reached for each 
wave function fmu and energy Emu- The total energy of the 
bound He state \miv\, 17121^2) can then be found, using either 
Eq. © or 



E — y ^ £mu — Edir — Ec 



(14) 
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m3;.3,m4fc(r-l,r2) = — [W^ma (^ix ) /mg ^3 (zi ) 



2.2 Continuum states of the helium atom 

The He state in which one electron occupies the bound 
(ma 1^3 ) orbital, and other occupies the continuum state 
(nuk) is denoted by \m3U3,m4k) = | Wma/mgi^a, Wm4/m4fc)- 
The corresponding two-electron wave function is 

1 

-W^m4(ri_L)/„4ft(2:i)W"m3(r2±)/m3i/3(22)] • (15) 

Here fm4k{z) is the longitudinal wave function of the con- 
tinuum electron, and k is the z-wavenumber of the electron 
at |z| — » oo (far away from the He nucleus). 

We can use Hartree-Fock theory to solve for the ion- 
ized He states as we did for the bound He states. Since the 
continuum electron wave function fmikiz) is non- localized 
in z, while the bound electron wave function fm^v^iz) is lo- 
calized around z = 0, it is a good approximation to neglect 
the continuum electron's influence on the bound electron. 
We therefore solve for the bound electron orbital using the 
equation 



2m„ dz'^ 



Ze Vmsiz) 



/"is "3 (z) 



-•'nsvsfmsiysiz) . (16) 



The continuum electron, however, is influenced by the bound 
electron, and its longitudinal wave function is determined 
from 



2me dz^ 



ZeWm^iz) 



+e 



dz \ f„i,_^i,^(z')\^ Dm^ruiiz — z') - Ef 



fm4k{z) 



^Z fm^kiz )fm^i^^{z )-Em^m4{z Z )fm^i^4{z) . 



(17) 

where Sf — Em^k = 7i^fc^/(2me). Here, the bound electron 
orbital ImaJ^s) satisfies the same boundary conditions as dis- 
cussed in Sect. l2.ll The shape of the free electron wave func- 
tion is determined by the energy of the incoming photon and 
the direction the electron is emitted from the ion. We will 
discuss this boundary condition in the next section. The to- 
tal energy of the ionized He state Imsz^a, mik) is simply 



E- 



(18) 



Note that the correction terms _Bdir and -Bcxc that appear 
in Eq. (|14p do not also appear in Eq. (|18p . The direct and 
exchange energies depend on the local overlap of the elec- 
tron wave functions, but the non-localized nature of the free 
electron ensures that these terms are zero for the continuum 
states. 



3 RADIATIVE TRANSITIONS 

We will be considering transitions of helium atoms from two 
initial states: the ground state, |00, 10), and the flrst excited 
state, 1 00, 20). 

In the approximation of an infinitely massive, pointlike 
nucleus, the Ha miltonian of the He atom in electromagnetic 
field is (see, e.g. jLandau fc Lifshit dTlQTTh 



j = l,2 



J = l,2 ^ 



.(19) 



where pj = —ihVj is the canonical momentum operator, 
acting on the jth electron, rj is the jth electron radius vec- 
tor, measured from the nucleus, and Atot{r) is the vector po- 
tential of the field. In our case, Atot{r) — Asir) + Acm{r), 
where j4s(r) and Acmir) are vector potentials of the sta- 
tionary magnetic field and electromagnetic wave, respec- 
tively. The interaction operator is Hint = H — Hq, where 
Ho is obtained from H by setting Aamir) — 0. The un- 
perturbed Hamiltonian Ho is responsible for the stationary 
states of He, discussed in Sect.[21 The vector potential and 
the wave functions may be subject to gauge transformations; 
the wave functions presented in Sect .[2] correspond to the 
cylindrical gauge AB{r) = x r. Neglecting non-linear 
(quadratic in Aem) term, we have 



Hin 



where 



2m„c 



[-n-J • Acmirj) + Aem{rj) ■ TTj], 



p + -Afl(r). 

c 



(20) 



(21) 



is the non-perturbed kinetic momentum operator: tt = 
rneT = me{i/h)[Hor — r Ho]- 

For a monochromatic wave of the form Aemir) oc ee"''', 
where e is the unit polarization vector, applying the Fermi's 
Golden Rule and assuming the transverse polarization (e ■ 
q = 0), one obtains the following general formula for the 
cross section of ab sorption of radiation from a given initial 
state \a) (see, e.g.. [Armstrong fc Nichol5ll972l) : 

a{Lj, e)=y—\e- {6|e'^"- j|a) |' 5{lj - tJta), (22) 

b 

where |6) is the final state, a; = gc is the photon frequency, 
^ba = {Eb — Ea)/h, and j is the electric current operator. 
In our case, j — (— e/me)(7ri -1- 772). 

We shall calculate the cross sections in the dipole ap- 
proximation - i.e., drop e"' '' from Eq. (|22|) . This approxima- 
tion is sufficiently accurate for calculation of t he total cross 
section as l ong as huj <^ nieC^ (cf., e.g., Potek hin fc Pavlov! 
1 19931 . Il997l for the case of H atom). In the dipole approxi- 
mation, Eq. (|22[) can be written as 



fbaS{0J - 



^ba 



where 



\{b\e ■ 7r|a)| 



2meiOl,a |,,| I >|2 

\{b\e ■ r\a)\ 



(23) 



(24) 



is the oscillator strength. In the second equality we have 
passed from the 'velocity form' to the 'length form' of the 
matrix element (cf., e.g., IChandrasekharl 1 1945^ . These rep- 
resentations are identical for the exact wave functions, but 
it is not so for approximate ones. In the adiabatic ap- 
proximation, the length represent ation [i.e., the right-hand 
side of Eq. lIMll] is preferable f see [Potekhin fc Pavlovlll993l : 



iPotekhin. Pavlov, fc Venturalll997l v 

To evaluate the matrix element, we decompose the unit 
polarization vector e into three cyclic components. 



e_e+ -I- e+e_ + eoeo, 



(25) 
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with eo = Bz along the external magnetic field direction 
(the z-axis), e± = (e^, ± \e.y)l\/2^ and = ■ e (with 
Q = ±,0). Then we can write the cross section as the sum 
of three components, 



(j(tJ, e) = (T+(ci;)|e+| + o-_(a;)|e_| + a()(Lj)|e()| 



(26) 



where a a has the same form as Eq. (|23|l . with the corre- 
sponding oscillator strength given by 



ha = T \Adba\ = \Mba\ , 

with 

Mba = {&|e* • r|a), 



(27) 



(28) 



where r = r / po and Uc = eB/{mec) is the electron cyclotron 
frequency. 

3.1 Bound-bound transitions 

Consider the electronic transition 

I a) = \mV,m2V2) = \ Wmfmv,Wm.^fm.^v2) 

— > |fe) = \m'u',m2V2) = Wm2ff„2„2). (29) 

The selection rules for allowed transitions and the related 
matrix elements are 



0-+ 



ao : Am = 0, Ai^ = odd, 
Am = 1, Az/ = even, 
Am = —1, A;/ — even, 

Affca = \f<Ti {gm' u'\frnij) {gm2V2\Sm2V-2) ^ 



fm2 U2 ) 5 



(30) 
(31) 
(32) 

where Am = m' — m, Av = v' — v. The oscillator strengths 
for bound-bound transitions from the states |00, 10) and 
1 00, 20) are given in Table [1] 

The selection rules (|30p - (|32p are exact in the dipole 
approximation. The selection rules in m follow from the con- 
servation of the 0-projection of total (for the photon and two 
electrons) angular momentum. Technically, in the adiabatic 
approximation, the y follow from the propert ies of the Lan- 
dau functions (e.g. JPotekhin fc Pavlovlll993l ). The selection 
rules in v follow from the fact that the functions 5^/^/ and 
jrav have the same parity for even v' — v and opposite parity 
for odd v' ~ V. 

In addition to these selection rules, there are approxi- 
mate selection rules which rely on the approximate orthog- 
onality of functions gm'v' and jrnv (for general v ^ v'\ Be- 
cause of this approximate orthogonality, which holds better 
the larger B is, we have 



(33) 



where |£| ^ 1 and e — > as Ap ±00. Therefore, the oscil- 
lator strengths for transitions with q = ± and Av = 2, 4, . . . 
are small compared to those with Ai/ — 0. The latter oscilla- 
tor strengths can be approximated, according to Eqs. (|27|) . 
dsn, (ESI and (O, by 



fba ~ 2(m + l)l^6aM, 

(a — Am — ±1, 1/' = u 



(34) 



smallness of matrix elements for transitions of the type 
\mv,m2V2) — > \m' v' ,m2v'2) with 1^2 7^ V2 for a = ± and the 
smallness of cross terms in the matrix elements of the form 
{gm2'^2\fmu) {gm' u'\fm2V2) wheu m' = m2 (i.e., the so-called 
"one-electron jump rule"); we have therefore excluded such 
terms from the selection rule equations above [Eqs. pop - 

mi 



3.2 Photoionization 

The bound-free absorption cross section for the transition 
from the bound state \b) to the continuum state ]/) is given 
by Eq. (|22p with obvious substitutions \a) \b}, \b) \f), 
and 



dfc, 



(35) 



where is the normalization length of the continuum elec- 
tron [j'^^^^^,^<iz\gmk{z)\'^ — 1] and k is the wave number of 
the outgoing electron (Sect. [221) • Therefore we have 



2ne^Lz 



{!(/. 



mech Ufbk 



k\e''"^e ■ Tr\b)\ 



TT\b)\ 



}■ 



(36) 



where k = y'^lm^ej /h and \f±k) represents the final state 
where the free electron has wave number ±fc (here and here- 
after we assume A; > 0). The asymptotic conditions for 
these outgoing free electrons are (cf., e.g., IPotekhin et al.l 
I1997I ) gm,k{z) ~ exp[i(/:)fc(z)] at z ^ ±00, where ipk{z) = 
\kz\ -\- (fcao)~^ In \kz\ and ao = /meC^ is the Bohr radius. 
Since we do not care about direction of the outgoing elec- 
tron, we can use for calculations a basis of symmetric and 
antisymmetric wave functions of the continuum - that is, 
in Eq. (f36)) we can replace {/fe| and by {/cvcn| and 

{/odd|. The symmetric state |/oven) is determined by the 
free electron boundary condition g'^f. even(O) = and the 
antisymmetric state |/odd) is determined by (?mfe,odd(0) = 0. 
Since the coefficients in Eq. (|17p are real, gmk,cvcn{z) and 

gmk , odd 

(2) can be chosen real. At 2 — » ±00, they behave as 
5mfc,{ovon,odd)(-2;) ~ sm[ip{z) -|- coustant] (where the value of 
constant depends on all quantum numbers, including k) . We 

A.\„ _ 



still have the normalization J^^' ^2 I5'"'''{°^'="'°dd) 



Similar to bound-bound transitions, we can decompose 
the bound-free cross section into three components, Eq. (|26p . 
Thus, using the dipole approximation and the length form 
of the matrix elements, as discussed above, we have for (q = 
±, 0)-components of the bound- free cross section 



0"bf,ci(ti') ~ — CTTh 



X |(/|e* • f\b) 



2£/ 
2 



pI 



(37) 



The same approximate orthogonality leads to the 



where |/) = |/ovon) or |/) = |/odd) depending on the parity 
of the initial state and according to the selection rules, and 
o"Th = (87r/3) (e^ /nieC^)^ is the Thomson cross section. The 
selection rules and related matrix elements for the bound- 
free transitions 

\b) = \miy,m2U2) = \ Wmfm,v ,Wm2fm,2V2) 
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> I/) = \m'k,m2V2) = \Wm'gm'k,Wm29m2i^2) (38) 

are similar to those for the bound-bound transitions [see 
Eqs. (O-dSil)]: 

(Jo : Am = 0, Au — odd, 

^2 '-'2 I /"^2 ) 1 (39) 

(j+ : Am = 1, Av = even, 

|/m^){S™/fcl/m2.V2)) , (40) 

cr_ : Am = —1, A;/ = even, 

M/f, = 0n [{g^'k\fm^){g 

~^m' !y,m2i^2 {9'm2^2 \ftrLu) {9 m' k\ fm2 1^2)) : 



(41) 



In this case, the condition Av = odd means that gm'fc and 
fmv must have opposite parity, and the condition Av — even 
means that gm'k and fmv must have the same parity. The 
oscillator strengths for bound-free transitions from the states 
|00, 10) and [00, 20) are given in Tabled 

Note that in Eqs. (|4Up and H41|) . the second term in 
the matrix element (of the form {5m2i/2 )(5™'fcl/m2i'2)) 
corresponds to transitions of both electrons. This appears 
to violate the "one-electron jump rule" and other approxi- 
mate selection rules discussed in Sect. 13. ll [see Eq. (|33|l ]. In 
fact, these approximate rules are not directly relevant for 
bound-free transitions, since the matrix elements involving 
a continuuum state are always small: (pm'fcj/mi/) ^ as the 
normalization length —* 00. Rather, we use a different set 
of selection rules to determine which of these 'small' matrix 
elements are smaller than the rest. The first is that 

{gra'k\fmv){g 

(42) 

when u'2 7^ U2. This selection rule is similar to the bound- 
bound transition case as {3m2i'2 l/™2f2) involves a bound 
electron transition, not a free electron transition. The second 
approximate selection rule that applies here is more compli- 
cated: terms of the form {gm'i,|/mi/)(5m2fc|/m2^2) are small, 
unless m' = m2 and U2 — v. This exception for m' = m2 
and V2 = V 11, due to the exchange term in the differential 
equation for the free electron wave function [Eq. p7|) ] , which 
strongly (anti)correlates the two final wave functions \gra'v) 
and |5m2fc)- If m-' = ni-2 and v = V2, then since {gm'u\fni2i^2} 
is not small (in fact, it is of order 1), {gm2k\f 71121^2) not 
be small but will be of the same order as other terms involv- 
ing the free electron wave function. In particular, the second 
selection rule means, e.g., that the matrix element for the 
transition from |00, 10) to jOO, Ok) is 

Moo,io^oo,ofc = (5ofc|,/io)(3oo|/oo) - (goo|,/io)(3ofc|/oo), (43) 

where the second term is non- negligible, but that the matrix 
element for the transition from jOO, 10) to |0fc,20), which is 



A^oo,io^ofc,2o ~ {g2o\.fio){gok\.foo), 



(44) 



is small compared to the other matrix elements and can be 
ignored (see Fig. [TJ . 

We make one final comment here about the effect of 
exchange interaction on the free electron state. If the ex- 
change term [the right-hand side of Eq. (fT7)) ] is neglected in 
the calculation of the free electron wave function, then the 
cross terms (i.e., those involving two-electron transitions) in 



the matrix elements of Eqs. (|40p and (|4H) are small and can 
be neglected. One then obtains approximate photoionization 
cross sections which are within a factor of two of the true val- 
ues in most cases and much better for ctq transitions. If the 
exchange term is included in Eq. (|17|) but the cross terms in 
the matrix elements are ignored, significant errors in the a± 
photoionization cross sections will result. To obtain reliable 
cross sections for all cases, both the exchange effect on the 
free electron and the contribution of two-electron transitions 
must be included. 



4 RESULTS 

Tables [T] and [5] give results for transitions of helium atoms 
from the ground state (|00, 10)) and the first excited state 
(1 00, 20)). Table [1] gives results (photon energies and oscilla- 
tor strengths) for all possible bound-bound transitions with 
Au ^ 1, for the field strengths B12 ~ 1, 5, 10, 50, 100, where 
B12 — 5/(10^"^ G). Transitions \a) \b) for a = — are not 
listed separately, being equivalent to transitions \b) — ^ \a) 
for a = +. One can check that the oscillator strengths fba 
presented in Table [T] for a = + are well described by the 
approximation (|34p . 

Table [2] gives results (threshold photon energies and 
cross section fitting formulas, see below) for all possible 
bound- free transitions. Figure [T] shows partial cross section 
curves for all bound-free transitions from the ground state 
of helium for B12 = 1. The transition [GO, 10) \0k,20) is 
an example of a 'weak' transition, whose oscillator strength 
is small because of the approximate orthogonality of one- 
electron wave functions, as discussed at the end of Sect. l3.ll 
It is included in this figure to confirm the accuracy of our 
assumption. Figures [2] and |3] show total cross section curves 
for a photon polarized along the magnetic field, for B12 = 1 
and 100 respectively. Figures 3] and [S] show total cross sec- 
tions for the circular polarizations, a = ±, for B12 — 1. 
Finally, Figs. [6] and [7] show total cross sections for a — i: 
and B12 = 100. 



4.1 Fitting Formula 

The high-energy cross section scaling relations from 
iPotekhin fc Pavlov! (|l993l ). which were derived for hydro- 
gen photoionization in strong magnetic fields, also hold for 
helium: 



Chi 



(1 



2mj+9/2 



2mi+7/2 



(45) 



(46) 



where m^ is the m value of the initial electron that tran- 
sitions to the free state. In addition, we use similar fitting 
formulae for our numerical cross sections: 

C 



O"bf,0 ' 



0"bf,± 



(l+^y)2-5(l-hS(yrTy- l))4(™.+i) 
, 0(1 + y) 

" (i + ^y)2-5(i + e(yrT^-i))4(-.+i) 



0"Th 



0"Th 



(47) 
(48) 



where y — Sf/huthi and hujthr is the threshold photon en- 
ergy for photoionization. These formulas have been fit to the 
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Figure 1. Partial cross sections ct^q^ _j versus final ionized electron energy for photoionization of the ground state helium atom 
((mi, 7712) = (1,0)). The field strength is 10^^ G. The transition |00, 10) — > |Ofc, 20) in the bottom left panel is an example of a 'weak' 
transition. We have ignored these transitions in our calculations of the total cross sections. 



cross section curves with respect to the free electron energy 
£/ in approximately the 1 - 10* eV range (the curves are fit 
up to 10^ eV for strong magnetic fields B12 = 50 — 100, 
in order to obtain the appropriate high-energy factor). The 
data points to be fit are weighted proportional to their cross 
section values plus a slight weight toward low-energy values, 
according to the formula (error in a) oc cre/"'^^. 

Results for the three fitting parameters, A, B, and C, 
are given in Table [2] for various partial cross sections over 
a range of magnetic field strengths. For photoionization in 
strong magnetic fields {B12 > 50) the cross section curves 
we generate for the cr+ and a- transitions have a slight 
deficiency at low electron energies, such that the curves 
peak at £/ ~ 10 eV, rather than at threshold as expected. 



These peaks do not represent a real effect, but rather re- 
fiect the limits on the accuracy of our code (the overlap 
of the wave function of the transitioning electron pre- and 
post-ionization is extremely small under these conditions). 
Because the cross section values are not correct at low ener- 
gies, our fits are not as accurate for these curves. In Tabled 
we have marked with a '*' those transitions which are most 
inaccurately fit by our fitting formula, determined by cross 
section curves with low-energy dips greater than 5% of the 
threshold cross section value. 
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Table 1. Bound-bound transitions \a) — » The photon energy hu!i,a = Ejj — Ea (in eV) and 
the oscillator strength for different polarization components a [see Eq. Il27l l1 . All transitions 
Ai/ ^ 1 from the initial states |00, 10) and [00, 20) are listed, for several magnetic field strengths 
Bi2 = -B/(10^^ G). The last two columns list the transition energies ?it<jj^ and oscillator strengths 
/^^, corrected for the finite mass of the nucleus, according to Sect. 15. Il 



Bl2 


cr 


\a) 


-.|6) 






fba 


ba 


f* 

Jba 


1 





|00, 10) 


|00 


11) 


147.5 


0.234 












-> 10 


01) 


271.8 


0.124 


— 


— 




+ 




^ 00 


20) 


43.11 


0.0147 


44.70 


0.0153 







|00, 20) 


|00 


21) 


104.4 


0.312 












20 


01) 


277.7 


0.115 








+ 




^ 00 


30) 


18.01 


0.00930 


19.60 


0.0101 








20 


10) 


100.7 


0.0170 


102.3 


0.0172 


5 





|00, 10) 


^ |00 


11) 


256.2 


0.127 




_ 








10 


01) 


444.8 


0.0603 








+ 




00 


20) 


66.95 


0.00459 


74.89 


0.00512 







|00, 20) 


|00 


21) 


189.2 


0.176 












20 


01) 


455.0 


0.0537 


- 


- 




+ 




00 


30) 


28.94 


0.00299 


36.88 


0.00381 








20 


10) 


151.1 


0.00512 


159.0 


0.00539 


10 





|00,10) 


|00 


11) 


318.9 


0.0974 


_ 


_ 








10 


01) 


540.8 


0.0457 


— 


— 




+ 




00 


20) 


79.54 


0.00273 


95.42 


0.00327 







|00, 20) 


|00 


21) 


239.4 


0.136 


_ 


_ 








20 


01) 


553.3 


0.0405 


_ 


_ 




+ 




00 


30) 


34.84 


0.00179 


50.72 


0.00261 








20 


10) 


177.0 


0.00301 


192.9 


0.00328 


50 





|00, 10) 


|00 


11) 


510.9 


0.0557 












10 


01) 


822.2 


0.0266 








+ 




00 


20) 


114.2 


7.85e-4 


193.6 


0.00133 







|00,20) 


|00 


21) 


396.7 


0.0776 












20 


01) 


841.1 


0.0235 








+ 




^ 00 


30) 


51.92 


5.37e-4 


131.3 


0.00136 








20 


10) 


246.5 


8.41e-4 


325.9 


0.00111 


100 





|00, 10) 


|00 


11) 


616.4 


0.0452 












10 


01) 


971.4 


0.0221 








+ 




00 


20) 


131.4 


4.52e-4 


290.2 


9.98e-4 







|00, 20) 


|00 


21) 


485.0 


0.0626 












^ 20 


01) 


993.4 


0.0195 








+ 




^ 00 


30) 


60.57 


3.13e-4 


219.4 


0.00114 








20 


10) 


280.7 


4.80e-4 


439.5 


7.51e-4 



5 FINITE NUCLEUS MASS EFFECTS 

So far we have used the infinite ion mass approximation. 
In this section we shall evaluate the validity range of this 
approximation and suggest possible corrections. 

It is convenient to use the coordinate system which con- 
tains the centre-of-mass coordinate item and the relative 
coordinates {vj} of the electrons with respect to the nu- 
cleus. Using a suitable canonical transformation, the Hamil- 
tonian H of an arbitrary atom or ion can be separated 
into three terms (IVincke fc BavdIlQSSI : iBave fc Vincyil990l : 
ISchmelcher fc CederbaumI Il99ll ): Hi which describes the 
motion of a free pseudo-particle with net charge Q and total 
mass M of the ion (atom), the coupling term H2 between 



the collective and internal motion, and Hz which describes 
the internal relative motion of the electrons and the nu- 
cleus. H\ and H2 are proportional to M~^, so they van- 
ish in the infinite mass approximation. It is important to 
note, however, that H^ (the only non-zero term in the in- 
finite mass approximation) also contains a term that de- 
pends on Mq^ , where Mo « M is the mass of the nucleus. 
Thus, there are two kinds of non-trivial finite-mass effects: 
the effects due to H\ -\- H2, which can be interpreted as 
caused by the electric field induced in the co-moving refer- 
ence frame, and the effects due to H3, which arise irrespec- 
tive of the atomic motion. Both kinds of ef fects have been 
included in calculations o nly for the H atom l|Potekhinlll994l : 
IPotekhin fc Pavlov! Il997l . and references therein) and He"*" 
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Table 2. Bound-free transitions |6) — > |/): The threshold photon energy /iLJthr (in eV) and 
the fitting parameters A, B, and C used in the cross section fitting formulas [Eq. Il48l l1 . All 
transitions from the initial states |00, 10) and |00, 20) are listed, for several magnetic field strengths 
Bi2 = B/(10i2 G). 



Bl2 


fj 


\b) 


- l/> 




rrii 




A 


B 


C 


1 





|00, 10) 


|00 


Ifc) 


1 


159.2 


0.96 


0.093 


1.43e6 








10 


Ofc) 





283.2 


0.89 


0.20 


8.83e5 




+ 




00 


2k) 


1 


159.2 


0.70 


0.061 


7.95e2 








10 


Ik) 





283.2 


0.86 


0.094 


1.30e3 




- 




-» 00 


Ok) 


1 


159.2 


0.62 


0.030 


8.89e2 







|00, 20) 


^ |00 


2k) 


2 


116.0 


1.00 


0.062 


1.78e6 








-> 20 


Ok) 





289.2 


0.88 


0.22 


8.71e5 




+ 




|00 


3k) 


2 


116.0 


0.66 


0.038 


3.94e2 








20 


Ik) 





289.2 


0.54 


0.14 


6.48e2 








|00 


Ik) 


2 


116.0 


0.62 


0.029 


5.82e2 


5 





|00, 10) 


^ |00 


Ik) 


1 


268.2 


0.86 


0.061 


8.39e5 








-* 10 


Ok) 





456.4 


0.69 


0.16 


4.60e5 




+ 




00 


2k) 


1 


268.2 


0.68 


0.036 


1.14e2 








10 


Ik) 





456.4 


0.83 


0.057 


1.93e2 




- 




-> 00 


Ok) 


1 


268.2 


0.60 


0.020 


1.36e2 







|00,20) 


|00 


2k) 


2 


201.2 


0.92 


0.039 


l.lle6 








20 


Ok) 





466.5 


0.65 


0.18 


4.39e5 




+ 




|00 


3k) 


2 


201.2 


0.65 


0.021 


5.95el 








|20 


Ik) 





466.5 


0.54 


0.084 


9.13el 








|00 


Ik) 


2 


201.2 


0.61 


0.015 


7.82el 


10 





|00, 10) 


|00 


Ik) 


1 


331.1 


0.82 


0.051 


6.58e5 








-> 10 


Ok) 





552.5 


0.63 


0.15 


3.51e5 




+ 




00 


2k) 


1 


331.1 


0.67 


0.029 


4.94el 








10 


Ik) 





552.5 


0.81 


0.046 


8.43el 




- 




^ 00 


Ok) 


1 


331.1 


0.59 


0.016 


6.00el 







|00, 20) 


^ |00 


2k) 


2 


251.6 


0.88 


0.033 


8.77e5 








20 


Ok) 





564.9 


0.59 


0.16 


3.31e5 




+ 




— > |00 


3k) 


2 


251.6 


0.64 


0.017 


2.64el 








|20 


Ik) 





564.9 


0.53 


0.069 


3.97el 








^ |00 


Ik) 


2 


251.6 


0.61 


0.012 


3.25el 


50 





|00, 10) 


|00 


Ik) 


1 


523.3 


0.73 


0.034 


3.74e5 








10 


Ok) 





834.2 


0.54 


0.11 


1.96e5 




+ 




^ 00 


2k) 


1 


523.3 


0.63 


0.020 


7.15e0 








10 


Ik) 





834.2 


0.77 


0.033 


1.22el 




- 




^ 00 


Ok) 


1 


523.3 


0.57 


0.012 


8.94e0 







Ino 20^ 


— > loo 


2k) 


2 


409.1 


0.79 


0.021 


5.02e5 








20 


Ok) 





853.0 


0.50 


0.13 


1.83e5 




+ 




^ 00 


3k) 


2 


409.1 


0.62 


0.0104 


4.04eO 








20 


Ik) 





853.0 


*0.52 


0.052 


5.88e0 








^ 00 


Ik) 


2 


409.1 


0.59 


0.0058 


4.13e0 


100 





|00, 10) 


^ |00 


Ik) 


1 


628.8 


0.69 


0.029 


2.96c5 








10 


Ok) 





983.4 


0.51 


0.101 


1.56e5 




+ 




^ 00 


2k) 


1 


628.8 


0.62 


0.019 


3.12e0 








10 


Ik) 





983.4 


0.75 


0.031 


5.33e0 








^ 00 


Ok) 


1 


628.8 


0.56 


0.012 


3.94e0 







|00, 20) 


-> |00 


2k) 


2 


498.0 


0.75 


0.018 


3.96e5 








20 


Ok) 





1008 


0.47 


0.12 


1.45e5 




+ 




^ 00 


3k) 


2 


498.0 


0.60 


0.0092 


1.81e0 








20 


Ik) 





1008 


*0.50 


0.050 


2.60c0 








^ 00 


Ik) 


2 


498.0 


0.58 


0.0042 


1.69e0 



10 Z. Medin, D. Lai, and A. Y. Potekhin 




Figure 2. Total cross section (tq versus photon energy for helium 
photoionization, from initial states {mi, 1112) = (1,0) (solid lines) 
and (2, 0) (dashed lines). The field strength is lO^^ G. The dotted 
lines extending from each cross section curve represent the eff'ect 
of magnetic broadening on these cross sections, as approximated 
in Eq. l|55ll, for T = 10*-^ K (steeper lines) and 10*^ K (flatter 
lines). 



Figure 4. Total cross section a+ versus photon energy for helium 
photoionization, from initial states {mi, 1712) = (1, 0) (solid lines) 
and (2, 0) (dashed lines). The field strength is lO^^ G. The dotted 
lines extending from each cross section curve represent the effect 
of magnetic broadening on these cross sections, as approximated 
in Eq. ([SSj, for T = lO" K. 
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Figure 3. Total cross section fro versus photon energy for helium 
photoionization, from initial states {mi,m2) = (1,0) (solid lines) 
and (2, 0) (dashed lines). The field strength is IQi"* G. The dotted 
lines extending from each cross section curve represent the effect 
of magnetic broadening on these cross sections, as approximated 
in Eq. l|55ll, for T = 10^-^ K (steeper lines) and 10*^ K (flatter 
lines). 



Figure 5. The same as in Fig. H but for ct_ 



ion ijBezchastnov et al.l[l99i : IPavlov fc Bezchastnovll2005l ). 
For the He atom, only the second kind of effects have been 
studied (|A1-Huiai fc Schmelchej|2003al lbh. 
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Figure 6. The same as in Fig.|4l but for B = 10^'* G. 



10 



XI 



0.1 



0.01 





"1 1 ' ' ' : 
B = 10i4 G : 








■■.\ '\ 




\ 

■■s \ — 
\ 

\ - 

\ \ 
'\ \ - 
'\ \ 

\ \ 

^ \ 
■\ \ - 

\ \ 



100 



1000 10000 

hcj (eV) 



Figure 7. The same as in Fig. \E\ but for (T_ 



5.1 Non- moving helium atom 

The state of motion of an atom can be described by pseu- 
dom omentum K , which is a conserved vector since Q = 
(e.g.. lVincke fc Bavdl 19881 : ISchmelcher fc Cederbaumlll99lh . 
Let us consider fi rst the non-moving helium atom : K = 0. 

According to lAl-Huiai fc Schmelcherl (|2003al ). there are 
trivial normal mass corrections, which consist in the ap- 
pearance of reduced masses me/(l ± rrie/Mo) in H3, and 
non-trivial specific mass corrections, which originate from 
the mass polarization operator. 



The normal mass corrections for the total energy E of 
the He state |mi!^i, 77121^2) can be described as follows: 



E{Mo,B) = 



E{oo,{l + m,/MofB) 
I + me /Mo 



3 



(49) 



where Q.^ = {me/Mo)uj^ (for He, hQ.c = 1.588Bi2 eV). The 
first term on the right-hand side describes the reduced mass 
transformation. The second term represents the energy shift 
due to conservation of the total z component of the angu- 
lar momentum. Because of this shift, the states with suffi- 
ciently large values of mi -I- m2 become unbound (autoion- 
izing, in analogy wit h the case of the H atom considered by 
IPotekhin et al.(ll997l '). This shift is also important for radia- 
tive transitions which change (mi 4- m2) by Am 7^ 0: the 
transition energy TiuJba is changed by 7if2cAm. The dipole 
matrix elements Mba are only slightly affected by the normal 
mass corrections, but the oscillator strengths are changed 
with changing u)ba according to Eq. (I27|) . The energy shift 
also leads to the splitting of the photoionization thresh- 
old by the same quantity 7if2cAm, with Am = 0, ±1 de- 
pending on the polarization (in the dipole approximation). 
Clearly, these corrections must be taken into account, unless 
Q.C <^ LUba or Am — 0, as illustrated in the last two columns 
of Table [U 

The specific mass corrections are more difficult to evalu- 
ate, but they can be n eglected in the considere d B ran ge. In- 
deed, calculations bv lAl-Huiai fc SchmelchMl (|2003al ) show 
that these corrections do not exceed 0.003 eV at B < 10'' Bo- 

5.2 Moving helium atom 

Eigenenergies and wave functions of a moving atom depend 
on its pseudomomentum K perpendicular to the magnetic 
field. This depen dence can be described by Ha miltonian 
components (e.g.. [Schmelcher fc Cederbaumlll99ll) 



EL 

2M 



^ Mc 

3 



K -(B xr^ 



(50) 



where is the sum over all electrons. The dependence 
on is trivial, but the dependence on the perpendicular 
component K± is not. The energies depend on the abso- 
lute value K^. For calculation of radiative transitions, it 
is important to take into account that the pseudomomen- 
tum of the atom in the initial and final state differ due to 
recoil: K' = K + fiq. Effectively th e recoil adds a term 
PC q into the interaction operator (cf. iPotekhin et al.|[l997l : 
iPotekhin fc Pavlov|[l997l ). The recoil should be neglected in 
the dipole approximation. 

The atomic energy E depends on K±^ differently for 
different quantum states of the atom. In a real neutron star 
atmosphere, one should integrate the binding energies and 
cross sections over the ifx-distribution of the atoms, in or- 
der to obtain the opacities Q Such integration leads to the 
specific magnetic broadening of spectral lines and ioniza- 
tion edges. Under the conditions typical for neutron star at- 
mospheres, the magnetic broadening turns out to be much 



^ For the hydro g en atom, this has been done by 
iPaylov fc PotekhinI lll995^ for bound-bound transitions and 
bv lPotekhin fc Pavlov! lll997l) for bound-free transitions. 
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larger t han the conventional Dop pler and coUisional broad- 
enings (jPavlov fc Potekhin|[l995l ) . 

At present the binding energies and cross sections of 
a moving helium atom have not been calculated. However, 
we can approximately estimate the magnetic broadening for 
T ^ |(A_E)min|/fcB, whcre {AE)min is the energy difference 
from a considered atomic level to the nearest level admixed 
by the perturbation due to atomic motion, and fce is the 
Boltzmann constant. In this case, the ifx-dependence of E 
can be approximated by the formula 



E{K^) = E{0) + 



Kl 
2M| 



(51) 



where E{0) is the energy in the infinite mass approxima- 
tion and M± = K±{dE/dK±)~^ is an effective 'transverse' 
mass, whose valu e iM± > M) depends on the quantum state 
consi dered (e.g., IVincke fc Bavd Il988l : iPavlov fc Meszaroi 
1 19931 ). 

Generally, at every value of K± one has a different 
cross section (j{uj, K±). Assuming the equilibrium (Maxwell- 
Boltzmann) distribution of atomic velocities, the K±- 
averaged cross section can be written as 



exp 



E{0) ~ E{K^ 



dE{K^ 



.(52) 



where -Emin ~ —fiLo. 

The transitions that were dipole-forbidden for an atom 
at rest due to the conservation of the total 2-projection of an- 
gular momentum become allowed for a moving atom. There- 
fore, the selection rule Am = a [Eqs. (|30p - (l32|l ] does not 
strictly hold, and we must write 



(53) 



where the sum of partial cross sections is over all final quan- 
tum numbers m' (with m' ^ and m' 7^ m2 for Az^ = 0) 
which are energetically allowed. For bound-bound transi- 
tions, this results in the splitting of an absorption line at a 
frequency ujba in a multiplet at frequencies uiba + Smflc + 



(M- 



-M-'^)Kl/2h, where Sm ; 



-m — a and M 1 



is the transverse mass of final states. For photoionization, we 
have the analogous splitting of the threshold. In particular, 
there appear bound-free transitions at frequencies uj < cjthr 
- they correspond to Sm < K'j_/{2M±hQc)- Here, ojthr is 
the threshold in the infinite ion mass approximation, and 
one should keep in mind that the considered perturbation 
theory is valid for K\/2Mi_ < |(Ai?)min| < TiLOtti^. Accord- 
ing to Eq. (|53p . a{u,K^^ is notched aX u < ojthr, with the 
cogs at partial thresholds uthr -\- 5mQ,c — K'j_/ {2Mj_h) (cf. 
Fig. 2 in IPotekhin fc Pavlov|[l997l ). 

Let us approximately evaluate the resulting envelope 
of the notched photoionization cross section (|53p . assum- 
ing that the 'longitudinal' matrix elements [{. . .) con- 
structions in Eqs. (|3Up - p2p ] do not depend on K±. The 
'transverse' matrix elements can be evaluated following 
IPotekhin fc Pavlovl IjlQQit ): in the perturbation approxima- 
tion, they are proportional to |^| ''''"'e"'^' where = 
A'ipg/(2ft2). Then 



a{Lo < uithr, K±) « a(ujthi, 0) exp 



M± LJthr — 



M flc 



2Mi 



(54) 



where 9{x) is the step function. A comparison of this approx- 
imation with numerical c alculations for the hydrogen atom 
(|Potekhin fc Pavlovl [l997i ') shows that it gives the correct 
qualitative behaviour of a{u!,K±_). For a quantitative agree- 
ment, one should multiply the exponential argument by a 
numerical factor ~ 0.5-2, depending on the state and po- 
larization. This numerical correction is likely due to the ne- 
glected ifx-dependence of the longitudinal matrix elements. 
We assume that this approximation can be used also for the 
helium atom. Using Eq. (|52p . we obtain 



a{Lu) ^ (T(tJthr) exp 



M Q.C 



(55) 



for u> < uJthv Here the transverse mass M± can be evaluated 
by treating the coupling Hamiltonian H2 a s a perturbation, 
as was done bv lPavlov fc Meszarod l|l993h for the H atom. 
Following this approach, retaining only the main perturba- 
tion terms according to the approximate orthogonality rela- 
tion (|33p and neglecting the difference between M and Mo, 
we obtain an estimate 



M 



E 



E 

6( Am — a) 



1 + LJba/^c 



(56) 



where \a) is the considered bound state (|00, 10) or |00, 20) 
for the examples in Figs. [2H7]) and \b} are the final bound 
states to which a = ± transitions \a) — > \b} are allowed. 
According to Eq. (|34p . the numerator in Eq. (|56p is close to 
m + 1 for a = + and to m for a — 

For the transitions from the ground state with po- 
larization a = — , which are strictly forbidden in the 
infinite ion mass approximation, using the same approx- 
imations as above we obtain the estimate o-_(aj) oc 
a+{iu)hD,ckBT/{kBT + hflcf. 

Examples of the photoionization envelope approxima- 
tion, as described in Eq. (|55p above, are shown in Figs.[2H3 
In Figs. |6]and[7] (for B = lO'^* G), in addition to the mag- 
netic broadening, we see a significant shift of the maximum, 
which originates from the last term in Eq. (|49p . Such shift 
is negligible in Figs. |4] and [5] because of the relatively small 
Qc value for B = lO" G. 

Finally, let us note that the Doppler and coUisional 
broadening of spectral feature s in a strong magne t ic fiel d 
can be estimated, follo wing JPavlov fc Meszarod l|l993l l. 
iPavlov fc PotekhinI (|l995l ) and lRaiagopal et al.l l|l997l) . The 
Doppler spectral broadening profile is 



with 



^/tFAcjd 



exp 



c V M 



2 ^'-^ -L 2 

cos ^fl-f-^^sin 9b 



-1/2 



(57) 



(58) 



where 9b is the angle between the wave vector and B. The 
coUisional broadening is given by 



A, 



1 



27r (lj - a;o)2 + (Acoii/2)2 



(59) 
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with 

?iAcoii = 4.8neaorcft ( ) 



1/6 



: 41.5- 



— T 



\ ao J 



(60) 



where Ue is the electron number density and rcff is an 
effective electron-atom interaction radius, which is about 
the quantum-mechanical size of the atom. The convolution 
of the Doppler, coUisional and magnetic broadening pro- 
files gives the total shape of the cross section. For bound- 
free transitions, the Doppler and coUisional factors can 
be neglected, but for the bound-bound transitions they 
give the corr ect blue wings of the spectral features (cf. 
IPavlov fc p"otekhin.,1995i ). 



6 CONCLUSION 

We have presented detailed numerical results and fitting for- 
mulae for the dominant radiative transitions (both bound- 
bound and bound-free) of He atoms in strong magnetic fields 
in the range of 10^^ — 10^'* G. These field strengths may 
be most appropriate for the identification of spectral lines 
observed in thermally emitting isolated neutron stars (see 
Sect.[l|. 

While most of our calculations are based on the infinite- 
nucleus-mass approximation, we have examined the effects 
of finite nucleus mass and atomic motion on the opacities. 
We found that for the field strengths considered in this paper 
{B < 10^"* G), these effects can be incorporated into the 
infinite-mass results to obtain acceptable He opacities for 
neutron star atmosphere modelling. For large field strengths, 
more accurate calculations of the energy levels and radiative 
transitions of a moving He atom will be needed in order to 
obtain reliable opacities. 
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